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Abstract. The source of cosmic gamma-ray bursts (hereafter GRBs) is usually be- 
lieved to be a stellar mass black hole accreting material from a thick disk. The mech- 
anism for the production of a relativistic wind by such a system is still unknown. We 
investigate here one proposal where the disk energy is extracted by a magnetic field 
amplified to very large values B ~ 10 15 G. Using some very simple assumptions we 
compute the mass loss rate along magnetic field lines and then estimate the Lorentz 
factor r at infinity. We find that Y can reach high values only if severe constraints on 
the field geometry and the conditions of energy injection are satisfied. We discuss the 
results in the context of different scenarios for GRBs. 

I INTRODUCTION 

Most of the sources which are now discussed to explain GRBs (the coalescence 
of two compact objects or the collapse of a massive star to a black hole (collapsar) 
[1-3]) lead to the same system : a stellar mass black hole surrounded by a thick 
debris torus. The release of energy by such a configuration can come from the 
accretion of disk material by the black hole or from the rotational energy of the 
black hole extracted by the Blandford-Znajek mechanism. The released energy 
is first injected into a relativistic wind and then converted into gamma-rays, via 
the formation of shocks probably within the wind itself [4,5]. The wind is finally 
decelerated by the external medium which leads to a shock responsible for the 
afterglow emission observed in the X-rays, optical and radio bands [6]. 

The production of the relativistic wind is a very complex question because of the 
very low baryonic load that has to be achieved in order to reach high values of the 
terminal Lorentz factor. Just a few ideas have been proposed and none appears 
to be fully conclusive. A first possibility to extract the energy from accretion is 
the annihilation of neutrino-antineutrino pairs emitted by the hot disk along the 
rotation axis of the system, which is a region strongly depleted in baryons due to 
centrifugal forces. The low efficiency of this process however requires high neutrino 
luminosities and therefore short accretion time scales [7]. Another possibility to 
extract the energy from accretion is to assume that the magnetic field in the disk is 



amplified by differential rotation to very large values (B ~ 10 15 G). A magnetically 
driven wind could then be emitted from the disk with a fraction of the Poynting 
flux being eventually transferred to matter. The energy can also be extracted from 
the rotational energy of the black hole by the Blandford-Znajek mechanism [8]. 
We present here an exploratory study of the case where a magnetically driven 
wind is emitted by the disk. Matter is heated at the basis of the wind (by vv 
annihilation, viscous dissipation, magnetic reconnection, etc.) and then escapes, 
guided along the magnetic field lines. Section II describes a "toy model" to explore 
the behavior of such a wind. Despite its extreme simplicity, we expect that it 
can help to identify the key parameters controlling the baryonic load. Our results 
are presented in section III and discussed in section IV in the context of different 
scenarios for GRBs. 



II A "TOY MODEL" 

We solve the wind equations with the following simplifications : (i) we assume 
a geometrically thin disk and a poloidal magnetic field with the most simple ge- 
ometry (straight lines making an angle 9 with the disk) ; (ii) we consider that 
a stationnary regime has been reached by the wind; (iii) we use non-relativistic 
equations (to obtain the mass loss rate we just need to solve them up to the sonic 
point, where v < 0.1c) but we adopt the Paczyhski-Wiita potential for the black 
hole 

ds GM bh 2GM BH 

$bh = with r s = 5 . (1) 
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We write the flow equations (continuity, Euler and energy equations) in a frame 
corotating with the foot of the field line, anchored at a radius r in the disk 

p v s(x) — rh , (2) 
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where x = £/r , £ being the distance along the magnetic field line, and p, P, e and 
v are the density, pressure, specific internal energy and velocity in the flow. 

The total acceleration g(x) includes both gravitational and centrifugal terms. 
In this exploratory study the power deposited per unit mass q(x) only takes into 
account the heating and cooling due to neutrinos. We assume that the inner part 
of the disk is optically thick (which is probably justified for compact object merg- 
ers but is more questionnable for collapsars except for low a-viscosity (a < 0.01) 
[10]). We include the following processes : neutrinos capture on free nucleons, neu- 
trino scattering on relativistic electrons and positrons and neutrino-antineutrino 



annihilation (heating); neutrino emission by nucleons and annihilation of electron- 
positrons pairs (cooling). The temperature distribution in the disk corresponds to 
a gemetrically thin, optically thick disk : 

/ r \3/4 l\- 

T v {r) = T* y—j ^- — ^/^ y ^* * s ^ e ^ em P era ^ ure a ^ r *) ' (^) 

The section of the wind s(x) is easily related to the field geometry because the field 
and stream lines are coincident. We adopt the equation of state computed by [9] 
which includes nucleons, relativistic electrons and positrons and photons. 

The acceleration g(x) along a field line is negative up to x = x\ for angles larger 
than 0i ~ 60° (60° is the exact value for a Newtonian instead of a Paczyhski- 
Wiita black hole potential). For x > x±, g(x) is dominated by the centrifugal force. 
The sonic point of the flow is located at a distance x s just below x\ (the relative 
difference never exceeds 1%). We solve the flow equations in a classical way by 
inward integration along the field line. We start at the sonic point by fixing trial 
values of the temperature T s and the density p s from which we get the velocity v s 
and the position x s (from the condition of regularity at x = x s ) and then the value 
of the mass loss rate m. We observe that at some position x cr , the velocity v begins 
to fall off rapidly while T reaches a maximum T max < T u (r ). We adjust T s and p s 
so that close as possible to and T max to T v (r ). 



Ill RESULTS 

We have studied the dependence of the mass loss rate rh on the different model 
parameters and found the following expression : 



m(r) ~ 3.8 10 



13 




r 

—Ar) 



g/cm 2 /s . (6) 



The geometrical function / is normelized in such a way that it is equal to unity 
for r = 4 r g and 6{r) = 85°. The very strong dependance of m with T u (r) (tenth 
power) is in agreement with what is found for neutrino driven winds in spherical 
geometry [11]. Figure 1 shows that m also strongly depends on the inclination 
angle. The other important parameters are the position in the disk and the mass 
of the black hole, while m depends only weakly on all other parameters like the 
size of the optically thick region (here r in = 3 r g and r out = 10 r g ). In the more 
general case where the source of heating is not restricted to neutrino processes but 
can also include viscous dissipation, magnetic reconnection, etc, we have obtained 
a very simple and general analytical approximation for m [12] 
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FIGURE 1. Mass loss rate from the disk for a constant (8 = 85°) and decreasing (from 90° to 
80° between 3 and 10 r g ) inclination of the field lines. The disk temperature is T* = 2 MeV at 
r* = 4 r g . The mass of the black hole is Mbh = 2.5 M©. 



where e is the rate of energy deposition (in erg/cm 2 /s) between the plane of the 
disk (x = 0) and the sonic point (x — x s ~ xi), A$ is the difference of potential 
(gravitational+centrifugal) between x = and x — x\ and 5 is a factor close to 
unity depending on the distribution of energy injection between x = and x = x s . 

We can now estimate the average Lorentz factor f = E/Mc 2 at infinity The 
total mass loss rate M and the power injected into the wind E are given by 

M = 2 / J M ,rir = 2.6 (^) (^) F g<0 g/s (8) 
and £ = 2 10 =. (5^) ( JL) ( 1Q51/ j; g/s/sr ) e Ig/s , (9) 

where F gco = J r ^° /'/ rg / [x; 0(x)] xdx is a function of the field geometry only; E 1 is the 
burst power in gamma-rays, Vt/An is the beaming factor and / 7 is the efficiency 
for the conversion of kinetic energy into gamma-rays. The wind is powered by 
accretion but at the same time the disk is heated by viscous dissipation and cools 
by emitting neutrinos. We assume that these losses represent a fraction a of the 
power E injected into the wind, so that we can estimate at r* = 4 r g : 
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From equations (8), (9) and (11), we can calculate the average Lorentz factor 

r = ( Mbh v ( ^ V s a2) 
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The value of F gco is 56 for a constant inclination 9 = 85° and 250 if 9 decreases 
from 90° to 80° between r = 3 and 10 r g . We therefore conclude that large terminal 
Lorentz factors can be reached only if several severe constraints are satisfied : (i) 
low -Fg eo values, i.e. quasi-vertical field lines; (ii) low a values, i.e. good efficiency 
for energy injection into the wind with little dissipation ; (iii) low value of Q/47T, 
i.e. necessity of beaming. With the more general equation (7) we can obtain 
another simple and useful constraint : if the power e deposited below the sonic 
point represents a fraction x of the total power e to t injected into the wind, we have 

mc 2 5x 

For r = 4 r g and 9 = 85°, we obtain x\ = 2.182 and A$/c 2 = 0.18 which implies 
that x should not exceed 10~ 3 to have T > 100 ! 



IV DISCUSSION 

This study is clearly limited by its crude assumptions. However the severe con- 
straints we get show how difficult it may be to produce a relativistic MHD wind 
from the disk. An optimistic view of our results would be to consider that this diffi- 
culty could just be a way to explain the apparent discrepancy between the observed 
rate of GRBs and the birthrate of sources in the collapsar scenario, most collapsars 
failing to give a GRB. A more pessimistic point of view would be to conclude that 
the baryonic load of such winds is never sufficiently low so that they remain non 
relativistic. If one choose to rely on the Blandford-Znajek mechanism to power the 
wind [8] it should however be checked that this process is not "contaminated" by 
frozen material carried along magnetic field lines coming from the disk and trapped 
by the black hole. 
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